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ABSTRACT 

C^- ! 

O ■ We present optical linear polarimetry in the line of sight to HH 135/HH 136. 

The polarimetry of the field stars reveals two populations: one corresponds to a 
foreground interstellar component; the other originates in the interstellar medium 
in the vicinity of the Herbig-Haro pair and, therefore, can be used to study the 
magnetic field in the star forming region. Its direction is aligned with the jet 
of HH 135/HH 136, which could be an indication that the interstellar magnetic 
! field is important in the outflow collimation. The interstellar magnetic field 

magnitude was estimated to be of order 90 jj,G. According to recent numerical 

Q\ ', simulations, an interstellar magnetic field of such strength can be important in the 

definition of the outflow direction. There is also evidence that the associated dark 
cloud has an elongation parallel to the magnetic field. Our image polarimetry of 
the extended emission associated with HH 135/HH 136 shows a centro-symmetric 
pattern pointing to the knot E of HH 136. Previous near infrared polarimetry 
traces a different illumination center, namely IRAS 11101— 5829 - the probable 
exciting source of the system. This discrepancy can be explained if the YSO 
emission is completely blocked in optical wavelengths and the dominant optical 
source in the region is the knot E, whose nature is uncertain. A discussion of the 

.^J ' spectral energy distributions of HH 136-E and IRAS 11101—5829 is presented. 

Subject headings: ISM: Herbig-Haro objects — ISM: magnetic fields — tech- 
niques: polarimetric — ISM: individual (HH 135, HH 136) 
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Introduction 



Magnetic fields are believed to play a crucial role in the physics of jets and outflows in 
young stellar objects (YSOs). The models presently proposed to expl ain the outflow en gine 
in low mass YSOs rely on a magneto-c entrifugally driven mechanism (jShang et al.ll2006l and 
references therein; iFerreira et al.ll2006l ). Whether the magnetic field also defines the launch- 
ing mechanism and properties of jets in high mass YSOs is still unclear. Some observational 
findings suggest that the formation of intermediate- to high-mass stars also proceeds via 
disk accretion as in its low-mass counterparts, powering similarly high ly collimated outflows 
rtMartf et al1ll993l : brooks et al1l2003l : bavis et all 120041 . iGredell 120061 ) . On the other hand, 
the interstellar (IS) magneti c field can be relevant i n the maintenance of jets, as is sug- 
gested by the simulations of iDe Colle fe Ragal (120051 ). From an observational perspective, 
Menard fe Duchend (120041 ). based on a sample of classical T Tauri stars, suggested that the 
objects with bright and/or long jets might have their disk axes parallel to the interstellar 
magnetic field. 



HH 135 and HH 136 are very luminous Herbig-Haro (HH) objects discovered by lOgura fe Walsh 



( 119921 ) . who presented optical imaging and spectroscopy of the sources. They are located in 



Easte rn Carina in the southwest portion of the dar k cloud DCld 290 .4+01.9 ([Hartley et al. 



19861 ) near the bright-rimmed HH region BBW 47 (IBrand et al.lll9861 ). The recently discov- 
ered infrared cluster #59 from iDutra et al.l (120031 ) is also coincident with the HH pair. A 
more complete description of the optical objects in this region is presented by lOgura fe Walsh 
( 119921 ). The estimated distance to the optical/infrared objects in this region is in the 2.7 - 
2.9 kpc rang e. DCld 290.4+0 1 9, wh ich has a size of 28' x 12', is included in the CO(J=1-0) 



catalogue of lOtrupcek et al.l (120001 ). This line has a well defined Gaussian shape with a 
FWHM of 1.0 km s" 1 and V LSR = -19.8km s _1 . 

Infrar ed observatio n were recently used to detect and study the physical properties of 
the H 2 jet (lGredelll2006l ). His H 2 and [Fe il] continuum-subtracted narrow band images nicely 
trace the line emiss io n morphology in the o utflow . A CO molecular outflow is also present 
(jOgura et al.l 119981 ) . IChrysostomou et al.l (120061 ) present imaging circular polarimetry of 
HH 135/HH 136 which suggests a helical magnetic field in the outflow. 

The emission kno ts of HH 135 and HH 136 are distributed in a practically straight line 
(jOgura fe Walshlll992l ). This could be interpreted as an evidence for two jets with a com- 
mon origin. However, both jets are dominated by blue-shifted components, which has led 



Ogura fe Walsh! (119921 ) to propose that each HH object had a different source. Subsequent 
infrared polarimetry of the extended emission associated with the HH obje cts has shown 



that they have a common illuminating source, namely IRAS 11101—5829 (ITamura et al. 



19971 ). The apparent contradiction of these two observations can be avoided by the scenario 
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proposed by lOgura et al.l (119981 ). In this picture, the HH136 jet is deflected by a molecu- 



lar cloud, changing from a red-shifted 



extremity (see Figure 5 of lOgura et al. 



et near the IRAS source to a blue-shifted one in its 



1998). 



IRAS 11101-5829 is a luminous (10^) YSO flOeura fc Walshl Il992h associated with 



molecular masers o f different spec i es (jBraz et al.lll989t iTe Lintel Hekkert fc Chapmanlll996 



Walsh et al.lll997l ). iTamura et al.l (119971 ) suggested that it is a Herbig Ae/Be stars encircled 



by a dust disk. In particular , the presence of a 6.7GHz methanol maser points to a high- 
mass YSO (IWalsh et al.lll997l ). The masers profiles have Vlsr i n the range -24 to -21 km s _1 , 
indicating a kinematic distance of approximately 3 kpc. This velocity is similar to that of the 
DCld 290.4+01.9, which suggests that the IRAS s ource and the da r k cloud are associated - 



Molec ular e mission in CO and C S are reported by IZinchenko et al.l (119951 ). iBronfman et al. 



( 119961 ). and lOgura et al.l (I1998I ). The millimetric contin uum image of this source shows 
evidence of more than one emission core (jHill et al.ll2005l ). T hese data also indi cate a total 
cloud mass of 230 M , consistent with the mass estimated by lOgura et al.l (119981 ) of 150 M 
using CO observations. 

In this work, we present a study of the magnetic field in the interstellar medium (ISM) 
around the pair HH135/HH136 using polarimetric optical data. Polarimetry of the opti- 
cal nebula associated with the HH135/HH136 is also obtained. A brief discussion of the 
IRAS 11101— 5829 and HH 136-6 sources is presented. In Section [21 we describe the polari- 
metric data and their reduction. The results and discussion are presented in Section [3J In 
the last section, we summarize our findings. 



Observations and data reduction 



The observations have been done in 2005 February 12 with the 0.60-m Boiler & Chivens 
telescope at the Observatorio do Pico dos Dias, Brazil, operated by the Laboratorio Nacional 



in 


~. — - — , 
Maealhaes et al. 


( 


1996|). r 


1973|; 


Magalhaes et al. 


1996) 



pixels. The above telescope and instrumentation gives a field-of-view of 10(5 x 10(5 (1 
pixel = / . / 62). The data have been taken with an Ra filter. Polarimetric standards stars 



(jSerkowski et al.lll975l ; iBastien et al.lll988t iTurnshek et al.lll990l ) have been observed in order 
to calibrate the system and estimate the instrumental polarization. The measured values of 
the unpolarized standard stars were consistent with zero within the errors. Measurements 
using a Glan filter were also performed to estimate the efficiency of the instrument. They 
indicate that no instrumental correction is needed. 
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The reduction has been done using using the IRAF0 facility. The images have been 
corrected for bias and flat-field. Counts in the ordinary and extraordinar y images of each 



object were used to calculate the polarization using the meth od described i n lMagalhaes et al. 



(1l984h . We have utilized the IRAF package PCCDPACK flPerevrall2000f ) in the polarimet- 
ric analysis. We obtained the polarimetry of around 1 600 objects in the field-of-view. The 
results are presented and discussed in the following section. 

The ordinary and extraordinary images of the extended emission associated with HH 135/HH 136 
did not overlap allowing image polarimetry to be performed. It has been done considering 
circular apertures of 2 pixels radius centered in points distant from each other by 4 

pixels (^2'.'5) in each CCD direction. The results are presented in Section [3j 

In addition, we have performed differential photometry using as calibrators USNO ob- 
jects in the image: they are 490 in total. With this we could estimate R magnitudes for all 
objects in the field. 



Results and discussion 



3.1. Magnetic field geometry 



The direction of the magnetic field component in the plane of the sky can be traced by 
the position angle of the optical polarization. It is valid if one assumes that the polarization 



originates from the dicroic absorption of the s tar light by non-spherica 



align ed by the superparamagnetic mechanism (IDavis fc Greenstein 



197ll ; a recent review on grain alignment can been found in 



1951 



Lazariar 



interstellar grains 



ianl 120031 1 



Purcell fc Spitzer 



Figure [T] shows the number distribution of the position angle of polarization, 6, for 
objects with P/a P > 5, which corresponds to cr e < 5?7. We have also discarded objects that 
have positions superposed to the outflow. Using these restrictions, we reduce our sample to 
303 objects. The distribution is clearly bimodal with peaks at approximately 55° and 100°. 
Therefore, we performed a two-Gaussian fit which is also shown in Figure [lj The fitted 
parameters and errors for bins of 10° are shown in Tabled] (lines 1 and 3). The results are 
statistically the same for smaller or larger bin widths. An inspection of the data shows that 
these two populations have distinct spatial distributions and polarization magnitudes. This 
is illustrated in Figures [2] and [3] in which we have respectively plotted the results for objects 



2 IRAF is distributed by National Optical Astronomy Observatories, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under contract with the National Science Foundation. 
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with polarization modulus smaller and larger than 1.5% (an arbitrarily chosen number). 
The objects with small values of polarization tend to be distributed in regions where the 
extinction is less pronounced (Figured- left). In the right panel of Figure [21 we show the 
histogram of the position angle of this sub-sample as well as a Gaussian curve with the 
same mean and dispersion of that of Figure [1] centered at 100°. The agreement indicates 
that the population responsible for this peak in Figure [1] is well represented by polarization 
magnitudes smaller than 1.5%. The large vectors tend to be located in a strip running from 
the southwest to the northeast of the image which roughly corresponds to the dark cloud 
(Figure [3]- left). Again, one of the Gaussian curves in Figured] fits well the distribution of 
position angles. 

Our interpretation of the above results is that the population with smaller values of po- 
larization corresponds to foreground objects in the line of sight to the HH pair, while the more 
polarized objects have their polarization produced by grains associated with DCld 290.4+01.9, 
hence tracing the magnetic field alignme nt in the star forming region itself. To test this hy- 



pothesis, we have used the compilation of lHeiled (120001 ) to verify the polarization behavior in 



a larger field-of-view. We selected the objects within a 5°x 5° field centered at HH 135 and 
with Pjo P > 3 (85 objects). The number distribution of the position angles and a Gaussian 
fit are shown in Figure HI The Gaussian parameters are presented in Table [1] (fourth line). 
In spite of the larger dispersion, the mean position angle of Heiles' objects compares well 
with that of our suggested foreground component. The mean polarization magnitude of 
Heiles' objects is 1%, which is also consistent with our data. These results corroborate the 
supposition that the population having a mean angle around 100° corresponds to the large 
scale, and probably foreground, ISM. 

Another way to constrain the origin of each population is to determine the behavior 
of the polarization with distance, which, however, cannot be properly estimated with our 
data. From a statistical point of view, a faint object is generally farther than a bright one. 
So, an alternative, despite limited, approach is to check the polarization dependence with 
magnitude. Figure [5] shows that the polarization tends to increase with magnitude. This 
corroborates our hypothesis that the larger polarization values are associated with objects 
at larger distances. 

The above discussion gives us confidence that the small polarization component is as- 
sociated with the foreground ISM in the direction of the HH pair. Consequently, we should 
subtract this component from the observations to obtain the IS polarization produced by 
aligned dust in the star forming region. To estimate a value to the foreground component 
we have averaged the polarization of the objects with the observed polarization smaller than 
1.5%. This totalizes 212 objects with a mean polarization of 0.59 ± 0.36% @ 93.4° (the 
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quoted error is the standard deviation of the distribution). This value was subtracted from 
our sample of 303 objects. (We would like to note that all the arithmetics has been done us- 
ing the Stokes parameters Q and U). The number distribution of 9 for objects with P/ap > 3 
is plotted in Figure El The parameters of the Gaussian fit is shown in the second line of 
Table [TJ This distribution, which should represent the magnetic field orientation in the 
HH 135/HH 136 region, is similar to the uncorrected, but not the same. The mean position 
angle is at 41?9 ± 1?2. 

The direction of the interstellar magnetic field found above can be compared with the 
geometry of the young stellar object, in particular, with the outflow direction. The jet 
position angle (from North to East - equatorial coordinates) has been estimated using the 
line joining IRAS 11101—5829 and a given line emitting knot. They are: HH 135, HH 136/A- 
B-D-H. There are other knots, but whose emission occurs mainly in the continuum, so they 
could not trace the jet. The adopted position angle for HH 136 is the average of its four 
knots. The resulting position angles are: 40?0 for HH 135; and 37?9 ± 0?2 for HH 136. 
Therefore, the component of the interstellar magnetic field in the plane of sky (^42°) and 
the YSO outflow are approximately parallel. 

An interstellar magnetic field aligned with the jet is the best c onfiguration for the prop - 



agation of the outflow in the ISM, as recently demonstrated by iDe Colle &: Ragal (120051 ). 
These authors conducted two-dimensional numerical simulations of clumps (which, in their 
models, represent time-dependent ejection from YSOs) propagating in a magnetized ISM. 
They found that jets moving parallel to the ambient magnetic field can propagate to much 
longer distances when compared with those that propagate per pendicular to the magnetic 



field. They claim that this could explain the correlation found by lMenard fc Duchend (120041 ) 
for classical T Tauri stars; namely, the bright and long jets tend to be parallel to the in- 
terstellar magnetic fields. The jet associated with HH 135/HH 136 has a projected size of 
approximately 0.5 pc and high luminosity, so in this object we could be seeing the e f fect o f 



a parallel IS magnetic field keeping the jet. On the other hand, IChrysostomou et al.l (120061 ) 
have found evidence of a helical magnetic field in the outflow of HH 135/HH 136 based on 
infrared circular polarization, which can concur to collimate the jet. The present evidences, 
however, cannot state unambiguously which magnetic configuration is predominantly acting 
as the main large-scale collimating mechanism in this high-mass YSO. 

T he emission li nes of H2 and [Fe 11] in HH 135/HH 136 indicate a fast, dissociative J-type 



shock (iGredell 120061 ) . It is evidenced by the different space distributions of these emissions. 
In a J-type shock, the transverse (relative to the propagation direction) magnetic field is 
small. So the magnetic field direction inferred from our large-scale measurements may be 
similar to that in the ISM in which the shock propagates. However, we should again recall 
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a possible helicoidal field in the outflow (jChrysostomou et al.ll2006l ). which would produce a 
C-type shock or a J-type shock with precursors. More observations in order to constrain the 
detailed shock physical conditions - as, for instance, the H 2 v=0 transitions - may be helpful 
in disantangling the magnetic field geometry in the outflow region. 

We could also ask if the geometry of DCld 290.4+01.9 has some correlation with the 
magnetic field. Figure d shows a 0?5 x 0?5 DSS2 Red image centered at DCld 290.4+01.9. 
The lines represent the contour plot of the flux at 100 /im from IRAS. HH135/HH136 can 
be seen in the lower right quadrant, northeast of the HII region BBW 47. The denser portion 
of the cloud - as illustrated by the obscuration at optical wavelengths and dust emission at 
infrared - seems to be elongated in the northeast-southwest direction. If this is true, the 
interstellar magnetic field, the HH outflow, and the cloud e longation are all nearly parallel. 
This configuration is similar to what occurs in Lynds 1641 (IVrba et al.lll988f ). 



3.2. Magnetic field strength 



The strength of the magnetic field in the plane of the sky, B, can be estimate usinj 



B = (47rp) 



1/2 



A9 



B 



(i; 



where: p is the mass density of the ISM; v, the one- dimensional turbulent velocity; and 
the dispersion of the m agnetic field direction. This expression was proposed by 
Chandrasekhar fc Fermil (119531 ) and relies on the equipartition of turbulent kinetic and mag- 
netic energies and isotrop y of the moti ons in the medium. The overall ideia behind this 
formula is still accepted (IHeitschl 120051 ) . notwithstanding different effects could lead the 
above equation not to be the best estimate of the actual field: large fluctuations of the 
magnetic field amplit ude; acting of the nonmagnetic forces on the gas; inhomogeneity of the 
interstellar material (jZweibell Il996l ) . Recent numerical simulations of polarimetric maps of 
molecular clouds indicate that this formula overestimates the magnetic field by a factor 2 



(lOstriker et al.ll200ll : IPadoan et al.ll200ll : iHeitsch et al.ll200li Irleitschll2005l : iMatsumoto et al. 
200d 



The value of A9b in the star formation region can be estimated by the standard devia- 
tion of the fitted Gaussian to the number distribution of the position angle of the intrinsic 
polarization (see second line of Tabled]). This number is, however, an overestimate of the 
dispersion of the magnetic field direction since it includes the observational error associate d 
with the 9 measurement. Following the procedure suggested by lPereyra fc Magalhaed (120051 ). 
we obtain a A9b value of 13?3. The turbulent velocity was considered that of the dark cloud 
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and measured by lOtrupcek et al.l (120001 ) as 1 km s . A total mass density of 1.4 x 10 g 
cm -3 has been estimated from the number density of H 2 presented in lZinchenko et al.l (119951 ). 
which was based on CS(J=2-1) measurements near the IRAS source. Considering a factor 
of 0.5 to equation (jTJ) - as discussed above - we obtain an interstellar magnetic field strength 
of 90 fiG. However, we would like to note that this value should be interpreted as the order 
of magnitude of the field. The reason is twofold. On one hand, the observational values 
used in the magnetic field calculation have their own uncertanties. On the other hand, the 
values of the mass density, magnetic field d ispersion, and turbulent velocity can be tracing 
different portions of the ISM. iHeitschl (120051 ) obtained that a single estimation of B with the 
above procedure can be in error by a factor of 7. In addition, we would like to note that the 
above estimate of B should be associated with the large scale ISM around HH 135/HH 136, 
not with the outflow region. This value is larger than that measured in the diffuse ISM of 



a few AtG, but it is in the range obtained for star forming regions (see, e.g.. iGonatas et al. 



19901 : IChrvsostomou et al.lll994h . 



Recently, iMatsumoto et al.l (120061 ) have studied the alignment of outflows with magnetic 
fields in clouds cores through numerical simulations. They found that the outflow tends to 
be aligned with the large-scale (> 5000 AU) magnetic field if the magnetic field in the core 
is larger than 80 /iG. Our above estimate of 90 fiG may be interpreted as the strength in 



the d ark cloud, i.e., the initial magnetic field before the collapse (the B„ of IMatsumoto et al. 



20061 ). So the enhanced magnetic field in the collapsing core that originated the YSO has 



probably exceeded the limiting value of 80 fiG, making plausible that in this region the 
geometry of the magnetic field in the original cloud determined the direction of the YSO 
outflow. 



3.3. Imaging linear polarimetry of HH 135/HH 136 

Figure [8] shows the imaging linear polarimetry of the emission nebula associated with 
HH 135/HH 136. The background image is from our data: the object is seen twice because of 
the beam splitting produced by the calcite block. Only measurements with Pjop > 10 are 
plotted. The vectors sizes show that the polarization can be as high as 30%. The position 
angles define a clear centro-symmetric pattern, typical of scattering, whose center has been 



calcu lated and coincides with knot HH 136-E (following the denomination of lOgura fc Walsh 



19921 ). From North to South, this knot is the second bright source in our image. The centers 



calculated using the data over the HH 135 region or the HH 136 region are the same. 



HH 136-E is the brightest point in the R band image, having a magnitude of 14.37 mag. 
The knot B, the second brightest, has a magnitude of 14.78 mag, which corresponds to a 
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flux 30% smaller than knot E. In both estimates we have used an aperture radius of 3 pixels 
(= 1?8). 

Our imaging R-band polarimetry indicates HH136-E as the illuminator center of the 
scattering pattern, so it is un equivocally associate d with the region. Previous K-band po- 
larimetry of the same region (ITamura et al.lll997l ) shows also a centro-symmetric pattern, 
but whose center is coincident with IRAS 11101—5829. The dominant source in the infrared 



region is NIRS 17 (ITamura et al.lll997l - see also Figure 5 of lGredelll2006l ). which is coincident 
with knot J. They suggested that the IRAS source is obscured from our view by an optically 
thick disk, which is evidenced by the "polarization disk", but illuminates the associated 
nebula through the optically thinner pole. 

The optical depth of a dusty medium grows from infrared to optical wavelengths. There- 
fore, in the R band, the disk around the IRAS source can be optically thick even at its pole, 
thus preventing any flux to escape. This could explain why IRAS 11101— 5829 is not the 
source of the optical light being scattered in the nebula. However, it remains as an open 
question the nature of the knot E. 

HH136-E is the brightest R-band source in the outflow region and is associated with 
the infr ared source NIRS 9 , whose infrared colors are consistent with a pre-main-sequence 
object (ITamura et al.lll997l ). It has a very strong o ptical and inf r ared continuum , bein g 
practically absent of [S n] , H 2 and [Fe n] emission (IGredell 1200a : lOgura fe Walsh! 119920 . 
This makes a Herbig-Haro nature quite improbable. In spite of the suggestion from Schmidt 
plates that knot E has an important Ha emission, no slit spectroscopy at its exact position 
has been done. The spectral energy distribution (SEP) of the knots HH 136-A , HH 136-B, 
HH 136-E, and HH 135 can be done using DENIS jThe Denis ConsortiuiiJl2005h and 2MASS 
(ICutri et al.ll2003l ) data, and our photometry. None of these sources are detected in longer 
wavelengths. The SED of knot E has a rising slope from I to K band. A black body fit to 
this curve provides a bolometric luminosity of L , which would correspond to a ZAMS 
star of ^2 M . 

To explain the non-trivial radial velocity structure of the emission knots, and considering 
a same excit i ng sou rce for HH 135 and HH 136 as indicated by previous K-band polarimetry, 
Ogura et al.l (119981 ) have proposed a scenario in which one of the jets coming out from the 
exciting source is deflected by a nearby molecular cloud. The region of zero velocity i s located 
around knots D, E, F, and G (Figure [H see also Figure 7 of lOgura fe Walshlll992l). In thi s 
region there is also a slightly increase of the 12 CO antenna temperature (lOgura et al.lll998l ). 
Besides, there is a MSX source between knot F and G, which could represent the point of 
collision. So another possibility to the nature of knot E would be the point where the jet 
collides with the molecular cloud. 



-10- 



A spectroscopic analysis of knot E as well as a detailed velocity study of the whole 
complex can probably shed some light in what is going on in this region and the true nature 
of this object. 



3.4. Spectral energy distribution of IRAS 11101-5829 



Figure [TU1 shows the SED of IRAS 11101-5829 based on literature data (see figure legend 
for the references). To estimate the bolometric luminosity of IRAS 11101—5829 we integrate 
a cubic spline to its SED which provides a value o f 1.32 x 10 L ffi a t a di stance of 2.7 kpc. 
This is in agreement with pr evious estimates from lOgura &: Walshl ( 119921 ) of 1.39 x 10 4 L Q 
and from IWalsh et al.l (119971 ) of 1.59 x 10 4 L . Both of them are based on IRAS data, but 
consider different corrections. The above luminosity can be used to const r ain th e stellar 
mass. Using the massive stars evolutionary tracks of iBernasconi fc Maederl (119961 ) for Z = 
0.02, we estimate an interval of 11 - 25 M for the mass of the central object. The higher 
masses are obtained if the object is very young, with a convective envelope. 

The SED of an embedded source contains more information than just the luminosity 
of the central object. It results from the reprocessing of the stellar flux in the circumstellar 
envi ronment. To e xploit this aspect, we have used the recently available grid of models 



of iRobitaille et al.l (I2006h to reproduce the observed SED of IRAS 11101-5829. We have 
concentrated on models whose parameters are: 



1.0 x 10 4 L Q < L < 1.4 x 10 4 L Q - see discussion above; 
11 M < M < 25 - see discussion above; 



*0 



• i = 81?4. Following lOgura et al.l (119981 ). we consider that the jets make an angle 
with the plane of sky of ~5° and that the di s k is p erpendicular to the jets. Among 
the inclinations provided by IRobitaille et al.l (120061 ). we chose this value as a good 
approximation to IRAS 11101—5829; 

• aperture = 100.000 AU. At a distance of 2.7 kpc, this corresponds to 37". This is the 
largest aperture provided by the models. We have used it in order to approximate the 
the angular resolution of IRAS of ~ 90". 



We have then 307 models which have been visually inspected. In doing this, we have 
selected the best 41 models for which we have calculated the x 2 values. We have found that 
the models #3020025 and #3007152 produced the smallest \ 2 ■ Table[2]shows the parameters 
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of the models - we ask the reader to see a complete description of them in iRobitaille et al. 



( 120061 ). This simple modelling provides a object with a mass of ~ 13 M , which puts the 



YSO near the ZAMS, with a age of 10 6 years. 

The above result is unexpected considering the presence of jets that are typical of a 
younger object. To discuss that issue, we would like to initially recall the optical knot J 
(=NIRS 17). It is located at 1"8 from the IRAS 11101-5829 center position. However the 
error ellipse of this source is 19 x 5", so it includes the optical/infrared source (see Figure 
IH]). In the IR and optical, the knot emission is dominated by a s trong continuum and does 



not have [S n] emission lines (lOgura fe Walsh! Il992l ; iGredell 120061 ). So a Herbig-Haro nature 



appears to be ruled out. On the other hand, the SED presents two maxima: in ^60/im and 
« 2/mi. The above modelling has implicitly assumed that its near IR portion is caused by the 
circumstellar disk emission in the observer direction and therefore it should have the same 
center position of the far infra r ed em ission. This might not be the case for HH 135/HH 136. 



As proposed by iTamura et al.l (119971 ). a possible geometry is one in which the near infrared 
YSO emission (produced in disk) is obscured from the observer direct view, but it can flow 
from the pole and illuminate the nebular material in the jet region. We propose that the 
near IR SEP (the knot J) is the YSO reflected light in the pole cavity, as seen in HH 46 
JPopitalE978l l In such if the YSO emission was isotropic, the knot J should trace 



the YSO's SED. However, we should recall that the YSO emits anisotropically. Supposing 
the adopted inclination is correct, the knot J should receive the emission from a smaller 
inclination which has a larger near IR component. So the SED for an inclination of ^81° as 
seen from a direct view, should have smaller fluxes at these wavelengths. This would result 
in models having properties of a more embedded, consequently younger, object. To do a 
proper modelling it would be necessary to know the 3D configuration of the knot J and the 
YSO. 



4. Conclusions 

We have presented optical linear polarimetry of HH 135/HH 136 and the nearby ISM. 
Our main results are listed below. 

• The polarization of stars mapping the magnetic field in the star forming region is 
nearly parallel to the Herbig Haro outflow. We suggested that the elongation of 
DCld 290.4+01.9 is also parallel to the magnetic field. 

• We estimated the interstellar magnetic field strength as of order 90 fiG. 
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• The interstellar magnetic field direction and magnitude are adequate to play a role in 
determining the outflow direction and keeping the jet collimation. However, a collima- 
tion by a helicoidal magnetic field in the jet region is not discarded. 

• The R-band nebula associated with the HH pair has a clear reflection pattern whose 
center is the knot HH136-E. It seems to be a star, but its nature could not be securely 
determined and deserves more observational effort. 

• A simple modelling of the IRAS 11101—5829 SED indicates a mass of 13 M and an 
age of 1 Myr, which is inconsistent with the presence of jets in the object. A less 
evolved stage might be found if the knot J is assumed to be produced by the reflection 
of YSO light in the surrounding material. 
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Fig. 1. — Histogram of the position angle of the observed polarization for fields stars with 
P/up > 5 in the line-of-sight to HH135/HH136. A two-Gaussian fit is shown, whose pa- 
rameters can be found in Table HJ 
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Fig. 2. — Polarimetry of field stars in the line of sight to HH 135/HH 136 with P/cxp > 5 
and P < 1.5%. (Left) The vectors represent the direction and magnitude of the polarization 
whose scale is presented in the upper right of the figure. The background image is from 
the DSS2/Red. The epoch of the coordinates is 2000.0. (Right) Number distribution of the 
position angle of polarization for the same sample. The full line is one of the Gaussian curves 
obtained in the 2-Gaussian fit of Figure [TJ 
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Fig. 3. — The same as Figure [2] for stars with P/op > 5 and P > 1.5%. 




Fig. 4. — Number distribution of position angle of polarization for objects in the iHeiles 
(|2000|)'s catalog within a 5°x 5° field-of-view centered at HH 135 and with P/crp > 3. A 
Gaussian fit is also shown - see Table [1] for the parameters. 
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Fig. 6. — Number distribution of position angle of intrinsic polarization for objects in the 
line-of-sight to HH 135/HH 136 with P/ap > 3 after the foreground polarization subtraction. 
A Gaussian fit is also shown - see Table [1] for the parameters. 



-22 - 



DEC 




llbilim llhl3ni llhlin llhLlm 



Fig. 7.— DSS2 Red view of a 0?5 x 0?5 region centered at DCld 290.4+01.9. The contour 
plot of 100 (Km IRAS is seen superposed. The epoch of the coordinates is 2000.0. 
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Fig. 8. — R-band imaging polarimetry of HH 135/HH 136. Only measurements with P/op > 
10 are shown. The polarization vectors, whose scale is presented in the top right corner, are 
superposed on our image. The two images correspond to the ordinary and extraordinary 
beams separated by the calcite block. The gray retangle marks the most likely position of 
the ilumination source. The coordinate scale is with respect to the right image and the 
vectors. The epoch is 2000.0. 
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Fig. 9. — Different wavelengths sources in the region of HH 135 /HH f 36. The contour plot 
is the [SII] image (continuum included) from iGredell (120061 ). We use different symbols to 
represent the data origin: see legend at left. 
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Fig. 10. — Spectral energy distributio n of IRAS 11101-5829. The observatio nal data are 



2MASS: Cutri et al. 



(2003 


). MSX: 


Egan et al. 


(2003) 



ISO spectrum from the ISO data archive ( http: //www.iso.vilspa.esa.es/ida/index.html ). The 
dotted and traced lines represent two models from iRobitaille et al.l (120061 ) (see text for de- 
tails). 
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Table 1: Parameters of the Gaussian fits to polarimetric data in the line-of-sight to 
HH135/HH136 



Suggested origin 


Mean 


Error 


Standard 


Error 


Data 




(°) 


(°) 


Deviation 

(°) 


(°) 




HH 135/HH 136 region - with foreground 


54.9 


1.4 


14.2 


1.0 


this work 


HH 135/HH 136 region - foreground subtracted 


41.9 


1.2 


13.8 


1.0 


this work 


foreground ISM - 10' x 10' field 


100.5 


0.4 


10.2 


0.3 


this work 


foreground ISM-5°x 5° field 


107.7 


2.3 


24.8 


2.9 


Heiles (2000) 
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Table 2: Parameters of iRobitaille et all fl2006[ ) models that best fit the IRAS 11101-5829 
SED 



Para xnpfpT T^p^rrintinn 


1 n p 

V ( 1 1 LLVj 




Model S 


30071 52 

WW 1 -L \J £-1 


3020025 


Stellar Mass (M^,) 


12.87 


13.01 


Stellar Radius (R0) 


4.47 


4.49 


Stellar Temperature (K) 


29,390 


29,545 


Envelope Accretion Rate (M /yr) 


1.04 x 10~ 3 


6.71 x 10~ 4 


Envelope Outer Radius (AU) 


1.00 x 10 5 


1.00 x 10 5 


Cavity Angle (degrees) 


31.8 


29.4 


Disk Mass (M ) 


3.64 x 10~ 2 


1.69 x 10~ 2 


Disk Outer Radius (AU) 


86.6 


159 


Disk Inner Radius (Rsub) 


1.00 


1.00 


Disk Inner Radius (AU) 


9.51 


9.53 


Scaleheight factor 


0.837 


0.996 


Disk Flaring Power 


1.073 


1.104 


Cavity density (cgs) 


1.25 x 10~ 20 


1.12 x 10- 20 


Ambient density (cgs) 


6.09 x 10- 21 


6.54 x 10~ 21 


Disk accretion rate (M Q /yr) 


5.72 x 10~ 6 


7.04 x 10~ 7 


Total A(V) along line-of-sight 


4.90 x 10 3 


0.99 x 10 3 


Total Luminosity (L Q ) 


1.38 x 10 4 


1.38 x 10 4 


Disk Scaleheight at 100AU 


6.84 


8.77 



